THE SHIMARINE FIRING 
PHASE IECISI*N 

JAMES FItANCIS TUCKER 





■ ‘ V. 



(l,ilifrruin 



» '• . T- 



i 

I 



'PenfUf bookbinding 

2216 CLEMENT STREET San Francisco 21, Calif. BA. 1^352 



DIRECTIONS 


FOR BINDING 


BUCKRAM 

COLOR NO. 


LETTERING ON BACK 
TO BE EXACTLY AS 




PRINTED HERE. 


FABRIKOID 

COLOR 


TUCKER 




1954 


LEATHER 

COLOR 

Letter in Gold 


THESIS 

T825 


OTHER INSTRUCTIONS 

Letter on the fro 


it cover: 


THE SUBlIARIf^E 


FIRING 


PHASE DECi; 


5 ION 


JAMES FRANCIS 

1 


TUCKER 



UNITED STATES 

NAVAL POSTGRADUATE SCHOOL 




SUBMARINE FIRING PHASE DECISION 

BY 

JAMES FRANCIS TUCKER 
CAPTAIN, UNITED STATES NAVY 




.-V. • 



I 




i 

i 

I 



I 

I 



Library 
U. S. Nava 



al Post4?raduatc 



Moritui ey. CaUfomia 

THE SUBMARINE FIRING 
PHASE DECISION 

by 

James Francis Tucker 
Captain, United States Navy 



Submitted in partial fulfillment 
of the requirements 
for the degree of 
MASTER OF SCIENCE 

United States Naval Postgraduate School 
Monterey, California 




1954 



Thesis 

-r-r , 



THE SUBMARINE FIRING 



PHASE DECISION 



9^c 



James Francis Tucker 



27379 



This work is accepted as fulfilling 
the thesis requirements for the degree of 



MASTER OF SCIENCE 



from the 

United States Naval Postgraduate School 



William Peyton Cunningham 
Academic Associate 
Operations Analysis Curriculum 



Approved: 



Academic Dean 




4 » 



f 

I 




I 

1 

I 

I 



i 




PREFACE 



As indicated in the title, this paper considers the problem of 
the Submarine Commanding Officer who must determine salvo size 
and spacing in making a torpedo attack on a target vessel. 

The paper presupposes a knowledge of certain concepts of 
the theory of probability. The terminology and notation used in 
discussing mathematical probability correspond to the expressions 
in general usage in the literature, such as, "The Theory of 
Probability", M. E. Munroe, 1951. The paper is offered parti- 
cularly as a consideration for possible classroom discussion by 
operations analysis students. 

The writer wishes to thank Professor Thomas E. Oberbeck 
of the U. S. Naval Postgraduate School for his original suggest- 
ions, assistance, encouragement and cooperation in the prepara- 
tion of this paper. 
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GLOSSARY OF NOTATION 



A a 
B a 
C a 
D a 

f.g 



h a 



1 a 



L a 



a. p. 



n a 
m 3 



R a 
S a 




N. 

J 



N 



o 




track angle 
target true bearing 
target true course 

domain or area of uncertainty associated with position of MOT 
a probability density functions 
hit probability 

aim point designator: i = 1,Z, ,n 

target length 
s. a aim point sequence 

number of aim points in an aim point sequence 
number of hits required for a specified degree of damage 
a number of T/O's fired at a single aim point. 

: minimum acceptable probability of a hit 
target range 
target speed 
time on station in days 
maximum salvo size 

! number of T/O's onboard on jth day of patrol; j a 0, 1,---,T 
= number of T/O's with which S/M begins patrol 
: expected frequence of attack on jth day of patrol 
angle on the bow (AOB) 

maximum expenditure of T/O's in one attack on jth day of patrol 
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CHAPTER 1 



ELEMENTS OF THE PROBLEM 
1, 1 Statement of the Problem 

Each time a submarine commander decides to make an attack 
he must make a complete firing phase decision. This consists, 
among other things, of deciding: 

( 1 ) 

(a) How many hits are wanted on the target. 

(b) The price in torpedoes he is willing to pay for these hits. 

Contained also in the firing phase decision is: 

( 2 ) 

(a) Selection of the final inputs into the problem. 

(b) Selection of a torpedo spread technique to insure hitting 
the target. 

(c) Control of the firing of the 'torpedoes (particularly as 
to time). 

This paper will try to give aid to (1) by an analysis of (2) 
based upon statistical inferences concerning 2(a) and the mathema- 
tics of probability concerning 2(b) and 2(c). 

The basis for decision will be derived in terms of probability 
of hits obtained for torpedoes eJ^pendsd. 

The handling of one attack is then to be correlated with the con- 
duct of a complete patrol based upon what is desired to be acobiq(plished 
by this patrol, as specified in the patrol orders. 
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1.2 Patrol Objectives 



A patrol is considered as taking station in a geographical area 
under one of two alternatives as to length of stay: 

1. No restrictions as to time. 

2. Stay on station a specified number of days. 

Each alternative is pursuant to patrol orders. The first envi- 
sages pursuing every opportunity to attack, and early expenditure 
of torpedoes relieves the submarine of further conunitments in the 
area. The second reflects a need to keep the submarine in the area 
for a definite additional reason, (e. g. reconnaissance, rescue, or 
to provide a constant menace to the enemy over a specified period 
of time, with availability of relief, say, a governing factor). 

It will be considered that submarine patrol orders contain the 
order to "Inflict Meiximum Damage on the Enemy"^ which may be 
interpreted to specify two clearly distinct tasks: 

A. "Insure" kill of each target attacked, 
or 

B. "Insure" damage to each target attacked. 

Here the word "insure" is used in a probabilistic sense, as are also 
the words "kill" and "damage". 

It is desired that the submarine commander do his best with his 
opportunities in any one of four categories in which his orders place 
him, namely Al, A2, Bl, B2, as shown in the following table. 



1. See Theodore Roscoe, Submarine Operations in World War I^, 1949: 
Sub Pac Operational Plan, p. 
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A B 



Kill each target attacked. 


Damage each target attacked 


^ Leave when torpedoes 


Leave when torpedoes 


expended. 


expended. 


Space the kills. 


Space the damages. 


^ Stay on station a specified 


Stay on station a specified 


length of time. 


length of time. 



In a type 2 patrol, i. e. . when A2 or B2 is assigned, note that it is 
desired that the commander spreadj^l his effort over the time speci- 
fied so that his ship presents, in effect, a continuous threat to the 
enemy. [This reasoning should also leave him with some torpedo 
capacity to protect himself in a tight spot throughout the duration 
of a type 2 patrol. ] 

It will be seen that in a type 2 patrol the submarine Command- 
ing Officer must have definite knowledge of the expected frequency 
of enemy target traffic through his patrol area upon which to base 
his estimate of expected frequency of attack. This implies that there 
will be sources of intelligence from which he can estimate this target 
frequency and also that he must have some measure of his own search 
ability. It will be assumed that he canv from previous sources of 
intelligence of whatever origin, (e. g. intelligence reports, previous 
contacts, previous attacks, etc.) make such an estimate of the ex- 
pected frequency of attack on any given day of his patrol. 
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1. 3 Firing Phase Objective 



For uniformity, we define the firing phase as the period from 
the time of the last observation (last set of target data inputs) until 
the time when the last torpedo of the salvo crosses the target track, 
real or extended. 

In any attack, the final firing set-up, determined by whatever 
means, is put into the torpedo data computer at a given instant. 

The input data to be scrutinized here are: True bearing, range, 
angle on the bow(or true course) and speed of the target. True 
refers to the compass rose only; the input data are to be consider- 
ed as a set of random variables in the probability sense. The 
possible errors of these input data are evaluated to predict the 
size and probability distribution of the location error of the target, 
i. e. the nature of an area (or domain) of uncertainty, say D , 
which is moving and expanding with elapsed time, and into which 
the torpedo spread is fired. In particular, it is desired to know 
how unit mass , representing the totality of probability of target 
position, described by the probability density function f(D) asso- 
ciated with D , is distributed over this area D at a certain 
elapsed time t after the final inputs are made. 

From knowledge of D and f(D) at any given elapsed time t , 
it is desired to predict the aim point sequence, spread spacing, and 
salvo size necessary to obtain; (1) at least one hit with a minimum 
.acceptable probability p^ when damage is wanted: or, (2) at least 
m hits, with a minimum acceptable probability , when m hits. 
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on the average, are needed to sink the target. 

The Commanding Officer then completes the firing phase deci- 
sion on the basis of the information cited above. If a type 2 patrol 
is being conducted, it is also necessary to specify the salvo size for 
a given attack that is compatible with the requirements of the length 
of stay on station. The following chapters discuss the elements of 
this decision, i. e. salvo size as a function of the choice of one of 
the four patrol objectives and "how good the set-up is". 
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CHAPTER 2 

ELEMENTS OF ONE-DIMENSIONAL CONTINUOUS 
PROBABILITY DISTRIBUTIONS 

2. 1 Scope 

Certain fundsunentals of probability will be discussed, using 
simple models, in order to develop some concepts, definitions 
and notation for later use. 

2. 2 Target Position 

Consider a target T of length ' L on course C^ . For conven- 
ience, C^ will be taken as the positive direction on the axis of a 
one -dimensional coordinate system. Suppose that, at the instant t, 
the MOT (middle of target) is known to be on the x-axis, but its 
exact position is uncertain. The position of the MOT , at the in- 
stant t , may be regarded as a stochastic variable with an asso- 
ciated probability density function* f(xjjt). The function f(xj;t) 
will be termed the target position density function (at time t ). 

In the examples considered below, it will be convenient to distin- 
guish between the cases where f(Xj;t) is uniform and where 
f(xj^;t) is normal. 

Hence, when f(xj;t) is uniform, set f(x^;t) = f^(xj;t) where 

( gL for ^ Exj i X i ^ Exj 

0 elsewhere 

and a is the value of x, denoting the mean of the distribution. 

O 1 ” 

(See Fig. 2a). In this context, target location error is taken to be Ex^. 
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When 

(2.2) 



f(xj:t) is normal, set £(^j;t) = 

1 






/2itI5~x 



Exp 



where 



X, -a 
1 o 



where - oo < Xj < «» 

iigain, a^ is the mean of the distribution but in this context the 
target location error is taken to be 25‘xj . 

2, 3 Torpedo Position 

Let x^ represent a point on the same x-axis where a hypothe- 
tical torpedo (hereafter called T/O) crosses the axis at the instant 
t . Then h(x 2 ;t) , the probability of a hit on the target may be 

computed for any point x^; i. e. 



(2.3) 



h(x2;t) 




f(xj;t) dxj . 



The function h(x 2 lt) is shown in Fig. 2 b for the case where the 
target position error EXj^ exceeds L. In this connection, set 
Exj = M*L where M is a constant, M >1. 

Suppose it is known that a T /O aimed at any point on the x-axis, 
say Xj will definitely cross the x-axis at the instant t , but the 
exact point of crossing is known only uncertainly as a stochastic vari- 
able with the associated conditional probability density function 

V I 

g(x2V,t Ixj). 

Then, P{xyt) , the conditional probability, of a hit on the target 
when the T /O is aimed at the p oint x^ is given by the formula 



1. Some attributes of the normal distribution are discussed in 
Appendix A 
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(2.4) P(x-yt) = 1 h(x 2 ;t) gix^;! jx^) dx^ . 

J-(70 

By way of illnstrating P(x^;t) , assume g(x 2 ;t|xj) is also a 
uniform probability density function such that: 

g(x2;t|x2) = 

In this context, the quantity Ex^ will be termed the T/O error. 
Assume also that f(xj;t) is given by (2. 1 ). Finally, let 
Exj ;> ^ Li , i. e. , ~ M*Li , Ex^ = N»L , with 

M > N 1 . The graph of P(x 2 »t} under the assumptions of 
this paragraph is shown in Fig. 2 c. 

Whenever the instant t is not essential to the discussion, it 
will be convenient to write 

f(Xjjt) = f(Xj) 

hfx^rt) = 

P(x2;t) = P(xj) 

etc. 

Since the primary interest of this chapter is the probability 
density function associated with target position, it will be assumed 
that the T/O error, Ex 2 , is zero. Under this assumption, 

(2.5) x^ = X 3 

|P(x 2 > = h(x^) 

2.4 Computation of hit probabilities. 

An hypothetical example is inserted at this point, primarily 
as a vehicle to explain the notation used in the paper. 




for 



2 ^^2 - ^2 ^3 



2 ^^2 



elsewhere 
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It is convenient to formulate first of all a system for designating 
the crossing points of T/O's on the ic-axis. By assumption [see (2. 5)] 
the crossing point of the T/O is the aim point of the T/O. The choice 
of aim points t and the probable consequences thereof is a basic 
objective of this paper. A definite choice of equally spaced aim 
points will be called an aim point sequence (a. p. s. ). Inasmuch as 
it will be assumed throughout this paper that all T/O's cross the 
axis at the same instant t , the a. p. s. is not concerned with the 
order of firing; it specifies only the number of aim points, where 
aimed, and the spacing between aim i>oints. Thus, let an a. p. s. 
be defined as follows: 

a^ = point marking the mean of f(x^) 

d ^1/2 distance between adjacent aim points along the axis 
n a^number of aim points in the a. p. s. 

See Fig. 2 d. 

Then, according as n is odd or even the aim points in an a. p. are: 





n odd 






n 


even 


^1 


= a 

o 




^1 


= 


a A 

o — ^ 


^2 


= r.a + 
o — 


2d 


^2 




a T d 
o 


^3 


o 

II 


2d 


.^3 


s 


a + 3d 

O “ 








^4 


s 


a T 3d 
o 



etc. 

Attention is turned forst to the case where a single T/O is 
placed at each aim point of an a. p. s. , with spacing 2d = L , 
a so-called simple salvo. Later, the case of placing more th|ua one 
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T /O at each aim point of’ an a. p. s. will be discussed. 

For the case of a simple salvo, let denote the probability 

of a hit associated with the i th aim point of an a. p. s. , i = 1 , . . . , n. 

n 

Let P,(n) = XT P. . 

^ i=l ^ 

i. e. , Pj(n) is the total probability of one hit for a simple salvo 
of n aim points. 

By way of illustration. Table ^.Igives values of P. and Pj^(n) 
for simple salvos in two cases. The first case, designated in the 
table by f^ , is for f(xj) uniform with Ex^ = 4L [cf. formula 
(2. 1)] and the second case, designated by f^ , is for f(xj) 
normal with = L [cf. formula (2.2)] 













Table 


I 








n = 


1 


n = 


2 


n = 


3 


n = 


4' 


a. p. 


f 


f 


f 


f 


f 


f 


f 


f 




u 


n 


u 


n 


u 


n 


u 


n 


^1 


. 250 


. 383 


. 250 


. 341 


. 250 


. 383 


. 250 


. 341 


^2 






. 250 


. 341 


. 250 


. 242 


. 250 


. 341 


^3 










. 250 


. 242 


. 250 


. 136 


^4 














. 250 


. 136 


^1 


. 250 


. 383 


. 500 


. 682 


. 750 


. 867. 


1. 000 


. 954 



In the case of simple salvos, with 2d = L, note that 
there are no gaps nor overlaps in the spread spacing. This 
insures fo-r any n , a maximum value of P^ in the case of sym- 
metrical target position density functions , such as and 
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Two possible criteria for the choice of n are as follows: 

(a) Let = minimum acceptable probability of a hit with 

one T /O. In this case, n = maximum i for which p^^ p^ , 

(b) Let minimum acceptable total probability of one hit 

for a. p. s. selected. In this case, n = minimum t for which 

P P 
*^ 1 — • 

By way of illustration, suppose p^ is fixed at 0. 2; then from 
Table I, for f^ say, it can be seen that, n = 4 would not be 
used, but instead n = 3 would be used. If, however, a high 

P was desired, say 0. 9 , then n = 4 would necessarily be 

o • ■ 

% 

used. The choice between (a) and (b) is a matter of preference. 
Hereafter, the choice of n will be given in terms of both 
criteria. 

The possibility of more than one hit with a salvo to T /O's 
arises if 

(1) more than one T /O is fire^ at each aim point of an 
a. p. s. with the distance between adjacent aim points 
exactly equal to target length, i. e. , 2d = L ; 
or (2) one T/O is placed at each aim point of an a. p. s. 

with the distance between adjacent aim points less 
than target length, i. e. , 2d < L : 
or (3) more than one T/O is fired at each aim point and 
2d is less than L . 

These possibilities will be considered in the order just enumer- 
ated and will be referred to as Cases 1, 2, and 3, respectively. 

I 
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These cases will be considered, under the simplifying assumptions 
that 

(I) each T /O crosses the x-axis exactly at its intended aim point. 

(II) all T/O's of a salvo may be considered to cross the x-axis 
at virtually the s$.me instant t , even if more than one 

T /O is fired at any aim point. 

In considering Cases 1, Z, and 3, let denote the 

probability of obtaining two hits with a salvo of size nxm in which 

m T/0*s are placed at each of the n aim points of an a. p. s. 

with spacing 2d between aim points. 

' Under assumptions (1) and (II) above* it follows that in Case 

1 there are m hits if there is one hit and hence that 

P^(n,m,L) = P,(n) ; that is» the value of P (n,m, L.) 
m i m 

is identical with the value of Pj(n) for a simple salvo. For 

example* referring to Table I, if 3 T/O's are aimed at a^ * 

then three hits would be obtained with probabilities . 250 and 

.383 * respectively* for the target position density functions 

f and f of Table I, Similarly* for a salvo of 6 T/O's* it is 
u n 

found that P£(3* 2* L) has the values , 750 and ,867 , respec- 
tively* for the target position density functions f^ and f^ of 
Table I. 

Turning to Case 2* the probability of more than one hit depends 
in part upon the amount of "overlap"^ between successive aim points 
of the a. p. s. To illustrate* consider a salvo of two T/O's fired 
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4 * 4 
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in an a. p. s. of two aim points with spacing 2d « 7 Then for a 
target position den'sity function with « L (as in Table I) 

Pj(2,l,-^L) = .547 
P2(2.1,^L) = .197 

Finally, in Case 3, it is clear that if m T/O's are fired at 
each aim point of the a. p. s. with 2d < L , there is a possibility 
of m , 2m, etc. hits. In this case it is possible to calculate 
Pj^(n,m, 2d) where j is a positive integer having values 1, 2, etc. 
up to some maximum value depending upon 2d. For a given salvo 
of nxm T/0*8, it is important to note that whenever n 1 and 
2d < Li , a penalty is paid, namely, reduction in the value of 
P^(n,m,2d) to achieve the probabilities Pj^(n,m,2d) where 
1 . 

The spacing 2d = L has been measured along the x-axis, 

coincident with the assumed target course C . Consider the 

^ o 

T/0*s on any track T other than normal to the x-axis. Liet 

o 

A denote the track angle between T^ and C^ (measured as in 
Fig, 2 e), where for simplicity adjacent tracks are considered 
parallel). To intercept an interval L on the x-axis, the tracks 
must be spaced L sin A apart. L sin A is called the effective 
target length. For the one-dimensional case, it is clear that the 
probability of a hit associated with a given aim point does not 
depend upon the track angle. To avoid consideration of target 
width, very sharp track angles will not be considered here 
(roughly, A < 15° or A < 165° is not considered). 
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f (x,;t) 
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1 /Ex, 






Figure 2a: fy(x^;t) 



I 










i 

|X. 1/2M 

1 \ 


t ^ 

1/M 

i y 1 
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IV 
1 ^ 

1 • ' 

1 1 1 






< — 


— J 

(M-l)L ' 

ML ^ 









Figure 2 b: h(x 2 ;t) where Ex^ = ML , M>1 



P(x2;t) 
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n = 1 ^ Aim Point 

n = 2 — X - - Sequence 

n=3 )!> ^ ^ (a.p...) 

” X ^ a X 

Figure 2 d: The a.-p. s. depends upon n , the number of aim points, 

and the probability density associated with target position. 



L sin A 




Figure 2 e: Aim point spacing along the normal to T/O track. 
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CHAPTER 3 



ELEMENTARY PROBABILITY DENSITY FUNCTIONS FOR A 
STOCHASTIC VARIABLE WHICH DEPENDS LINEARLY ON TIME. 

Suppose that at a certain instant of time t , measurements 
are made of the position and speed of a point which is moving 
along a straight line at constant speed. Let the position and speed 
measurements both be subject to error and suppose that these mea 
surements may be regarded as observations of two independent 
stochastic variables. Then it is frequently necessary to "predict" 
the position of the point at some later time t + T. In this con- 
text, "predict" means to determine the probability density func- 
tion associated with the position of the point at time t + T. 

The future position of the point is a stochastic variable which 
depends linearly on time. If X denotes the future position of the 
point at time t + T and if x denotes the position of the point at" 
time t , then 

X = X + s 'T 
where s is the speed. 

This chapter discusses the probability density function asso- 
ciated with X when the probability density functions of x and 
8 are (1) both uniform and (2) both normal. 

3. 1 Probability Density Functions of x and s Both Uniform. 
Given ; (a) A uniform probability density fj * fj(x;t) in 

the interval (x^ - Ex * ^ 7 associated with the 
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measurement of the position of the point on the line at time t . 

(b) A uniform probability density f^ ^ in the 

interval (s - , b -i-i-E) associated with the speed of 

the point on the line* When the position of the point at time t is 
known exactly and when time is considered as being measured 
without error, then the position of the point after an elapsed 
time "Y is uncertain and the associated probability density 
function, here denoted by f^ » f 2 (sT;t) , is uniform in the 
interval (s^T' - ®o^ ^ 

Clearly at T' ^ 0 i f ^(x;t) completely describes the pro- 
bability density associated with the point, with mean at x^ 
and density f ^ in the interval (x^ - Ex, Ex) 

and 0 elsewhere. This would be the case at any instant t 
when a measurement of x is made. If no furl^her measurement 
of X is available, but the position of the point at time 
tj = t + ^ is required, then it is necessary to consider the new- 
stochastic variable 

X = X + s “t' 

The probability density function f^ of the new stochastic variable 

1 

is obtained by the convolution of and . 



Thus , 



= f3(X;tj) = fj(X;t) ♦ f2<X;t) 



= " ®T) f2(sT)rds 



For the simple case here, the distribution diagram for f^ is 
represented in Fig. 3 a. 

1. Cf: M. E. Munroe', Theory of Probability , 1951, Section 38. 
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It is geometrically evident that: 

(1) The mean of X is moving, and at time tj is at 

(2) The values of X are restricted to lie within the interval 

t*o + ®o”^ “ 4®* ^ -2®8T), + s^r + (^Ex + ^E*r)J 

which lengthens withT^ , 

(3) The area of the trapezoid representing the probability density 
function f^ over this interval must always equal one (the total'^ 
ity of probability). Consequently, the ordinate values, which 
represent probability density, must necessarily be reduced as 

^ increases. 

3. 2 Probability Density Functions of x and s Both Normal . 

Consider now. the case where each of the density functions f^ , 
^2 * ^2 uniform, but normal. In accordance with the nota- 

tion in Appendix A 

fj s ^i(x:x^, ^;t) 

Then for any time tj = t + T, the stochastic variable 
X s: X + s T has a normal probability density function which, 
by the addition theorem for normal distributions may be repre- 
sented in the notation of Appendix A 

Where X = x + s 

o o o' 

X X s 

The consequent movement of the normal variate X , the 
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expansion of and the thinning out of f^ is also clear in this 

case as shown in Fig. 3 b. 
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CHAPTER 4 

THE AREA OF UNCERTAINTY OF TAJp.GET POSITION AND ITS 
ROLE IN THE SUBMARINE FIRING PHASE DECISION. 



4, 1 Prelimliutry Considerationt and Notation 

Let the independent stochastic variables (R»B) describe the 
position of a point in a plane; e.g . , an observation of the range R 
and the true bearing B of a target. Let the uncertainty of the 
exact value of R (at time t ) be associated with the uniform 
probability density function 



Let B be associated with the uniform probability density 
function 



For definiteness, consider B measured in radians, say from a 
stationary submarine whose position will be considered a fixed 
datum point of the problem in any coordinate system desired. 

The joint probability density function associated with a pair of 
values (R,B) at t is then 



Thus, fj is positive over the annular sector shown in Fig, 4 a 
as the domain D^ s D^(t) . 




for R 



elsewhere 



b 





elsewhere 



for B^ - i-EB^ B i B^ + i EB 

O M O M 



f3(R,B;t) « fj(R;t) f2(B;t) 



Suppose now that no rrlatter what' point in is selected, there 
is a vector uncertain as to direction and speed, that can send 

this point (R,B;t) into a new position, say {P,jS^t$T) after an 
elapsed time 'Y « As heretofore, consider 'T as stop-watch time 
measured from the instant t without error. As shown in Fig. 4 b 
consider the angle that makes with B to be an indepen- 
dent stochastic variable, ^ith associated probability density 
function 



f^(e;t) s 




for- e 



1 



E6> i e s + 



E6 



O 7 

0 elsewhere 

Consider the speed S , i. e. , the magnitude of , to be an in- 
dependent stochastic variable, with an associated probability 
density function 



f5<S;t) 



for S 



jES < S 5 ^ 



E? 



0 elsewhere 

Then at the end of a specified time ^ , the probability density 

associated with the magnitude of the vector due to the variate 

S is, say /^1 o 1 ^ ^ c .1 

f5(ST:t) = 







BST 



0 elsewhere 

Evidently, for each possible origin (R, B;t) in , the vector of 
uncertainty has a domain of termination, say 

» D^(R, B;"*^, Y It) consisting of the annular sector 
shown in Fig. 4 b. And with is an associated probability density 
function 

f^(S,0;r;t) = f^(e;t) fgCSnt) 
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But V'T may originate from any point in and further may 

assume any one of the spt of values allowed by . It may be cbnsi 
dered that "maps" domain * D^(P,^ :t,T) 

where P and are new stochastic variables. Functionally, 
may be represented as 

(4 1) D 7 = [R.fj{R);B,f2(B);e.f4(0);S,f5(S):r;t] 

The region may be obtained as the result of a simple geome- 
trical construction and ig an eight-sided figure (see Fig, 4 c). 

The boundary of the area is comprised of the four sides of 
and the four sides of . With is associated a positive 
valued probability density function f^ = fy(P,^;t,7*) . The 
mathematical problem is to determine exactly how unit mass, 
representing the totality of probability, is spread over Dy*. fy 
denotes the description of this mass distribution over . 

No attempt to develop the analytic theory for the evalua- 
tion of the joint density function will be made here. Clearly, the 
probability density f^ of a set of values of the new variates 

at time (t+T) is zero at the boundary edges of D,^ and 
rises continuously, with possibly a different slope from each 

edge to a plateau of maximum probability density, say 
* 

fy S max fy 

It is evident that there exists no simple way of describing the 
infinite variety of configurations of and the associated proba- 
bility density function iy for and iy depend upon the rela- 

tive magnitudes of all of the parameters indicated in formula (4. 1). 
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For example, Fig. 4 d shows the regipn which results when , 

the area of uncertainty of target positipn at time t , is propogated 

with no speed error and a discrete course error. 

As a point pf reference in , define fo be the center 

of gravity of the mass distribution represented by f^ . The coordln- • 

ate s of A^ will be designated when the 

time parameters are essential. When P and ^ are independent 

stochastic variables then ^ = f* and ~ ^ where ^ and ^ 

denote mean values of P and ^ respectively. 

Alternatively, if f^ has a maximum value at a unique point 

A^ of Dy , then this point may be selected as a point of refer- - 

ence of O., and in this case the coordinates of A will be 
I o 

designated (r , p ) where 



In the familiar special case where fj is known (or approximated) 
in the form 



(4.2) £*(x,y) = 






Exp. 



/' (x-x 
' o' 



2<r^ 



ly-y^) 
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2 \ 



/ 



then A = (See Appendix A, Fig. A-2). In this special 

case the iso-density contours defined by = t where ^ is a 
constant are concentric ellipses. 

Proceeding frcm the formulation of and f^ with the 

simplest of assumptions concerning the component stochastic 
variates R,B,&,S, it would be desirable to turn to the case 
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where the component errors are independently and normally distri- 
buted. The reader will recognize the inherent matkcmatical difff- 
culties which are beyond the scope of thi#' paper. 

4. 2 The Submarine Fire Control Problem 

The concepts and notation of Section 4. 1 will now be related to 
the submarine fire control problem. Thus, let "(R,B;t) denote 
the actual range and bearing of MOT at the time of final observa-' 
tion of the target (instant of last inputs into Position Keeper of 
TDC before firing). and are the mean or expected ^|alues 

of R and .B at time t with respect to the distributions f^(R;t) 
and f 2 (B;t) , respectively. The functions fj and f^ are 
assumed to be known. It will be assumed further that the values 

of R and B are identical with the final observations of R 
o o 

and B (inputs into TDC). Similarly, B and S denote the AOB 

(angle on the bow) and the speed of the target at time t . 6^ and 

S^ are the mean or expected values of 0- and S at time t with 

respect to the functions f^(-0;t) and f 5 (S;t) respectively. Again, 

the functions f^ and f^ are assumed to be known and the values 

of 6 and S are assumed to be identical with the final obser- 
o o 

vations of B and S (inputs into TDC). Then the region D^ 
may be called the area of uncertainty of MOT at the time of 
last observation of the target. The region D^ then represents' 
where the MOT may be found at time t+T (yunits of time after 
the last observation) when there are uncertainties in target 
range, bearing, course, and speed. The region D^ will be 
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termed the area, of uncertainty of MOT at t+ 7^ . 

The intervj^l may be considered as the sum of two intervals 
and Xz • 

The interval represents the "dead-time" between the last 

observation of the target and the instant of firing: the interval 

• which is a function of the time required for a "perfect" 

T/O to reach when the T/O is fired at the 

instant t+Tj • Briefly, Xz T/O running time and 

Tj = T2(T,) 

The submarine is assumed to have a TDC which can predict 
within the interval X\ » 1* e. , T/O running time is calculated 
before the T/O must be fired. 

To this point, the direction of the target's motion has been 
discussed in terms of 0 , the AOB. It will be convenient to 
speak now of the target's course C , where, of course, 

C = C(0,B) and in particular -■ C(0^,B^). Inasmuch as 
is the target's course as it appears in the TDC in accordance 
with previous assumptions, it is convenient to consider the aim 
points of the aim point sequence (see Ch, 2) of the T/O salvo to 
be along the course line extending through the point a^ . 

Unlike the simple one -dimensional case discussed in Chap- 
ter 2, there is in the two-dimensional problem no simple way 
to calculate the probability of a hit P. associated with the ith 
aim point of the a. p. s. when f^ is generated by f^, f 2 » f^, 
and fg . To discuss the two-dimensional case further, it is 
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convenient to use the simplest of continuous probability density 



= f(x,y) = 

With this simplification, is a rectangle which will be 

designated D(t+7^) . See Fig. 4 e. 

Recall that D represents the area of uncertainty of position 
of MOT at time t+ 7^. The tinie interval is computed for the 

collision instant of a "perfect” T/O with the point a . Let C • 

- ' ' ■ o o 

the "most likely" target course, through a^ ,be identical with the 
x-axis of coordinates. Let the collision instant t+7^ = t^ . To 
simplify further, assume that all T/O's have infinite speed through 
D(to) , and that all T/O's of a salvo pass through D(t^) at the 
same instant t^ . The width of the target will be ignored. For 
T/O's traversing paths parallel to the y-axis, the probability for 
a hit at a given aim point along would be 



Referring to Fig. 4 f, it will be seen that t^e probability of a hit 
for such a track is simply the volume under f(x,y) "swept out" 
by the T/O "path" of width L sin A . 



functions as a model for f.^ . 



To simplify them, consider 





which reduces to the result (2. 2) for the one -dimensional case of 
Chapter 2. For track angles other than 90° this is not the case. 
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Consider now the assumption that the T/O has infinite speed in 
the light of the simple example just discussed. Taking into account 
the finite speed of the T/O, a very precise calculation of probabi- 
lity of a hit would require determination of the earliest possible 
instant, say alter firing at which the T/O enters a 

region D(t^- Sj) in which f(t^- 5j) has a positive value and the 
latest i>ossible instant, say at which the T/O leaves a 

region D(t^+S^) in which f(t^+S‘ 2 ). has a positive value. When • 
f(t^- 6 "j) and f(tQ+S' 2 ) differ sufficiently over the "path swept 
out" by the T/O, there may be an appreciable error in assum- 
ing that the T/O passes with infinite speed through a region 
D(t^) and sweeps out a path under, the probability density func- 
tion f(t^) . Moreover, the statements above that the T/O 
enters D(t^-S^j) and leaves imply that these regions 

have certain boundaries which will not be true when f has 
positive values over the entire plane as, for example, when at 
least one of the component probability density functions f j, f 2 » 
f^ , or fg is normal. Under these circumstances it may be 
convenient to consider a region, say D(t^-5'^) bounded by an 
C-valued iso-density contour 'of f(t^- (cf. remarks follow- 
ing formula (3. 2)) which the T/O enters at another 

region, say D(t + ?,) , bounded by an ^-valued iso-density 

contour of f(t +5".} which the T/O leaves at t +5*, and deter- 
o d o ^ 

mine whether f(t -S~,) and f(t +^ 5 ) differ "significantly" 

Ox O M 

over the path swept out by the T/O. When these two functions 
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# m 



do not differ significantly, th«n the probability of a hit associated 
with an aim point may be visualized as the ojolume under f(t^;5) 
swept out by the path of the T/O through D(t^;£) . 

In actual practice, the probability density functions f^, • 

f^ , and fg associated with target range, bearing, course, and 
speed inputs may not be known. However, it seems evident from 
the analysis of this paper that these functions are essential opera- 
tional data which are necessary for the submarine fire control 
problem if there is to be a rational basis for T/O salvos. It is 
likely that all four functions are normal density functions with 
standard deviations which maybe estimated by operational trials. 

Appendix B gives a brief discussion of the numerical approxi- 
mation of by high speed digital computers when • ^2 * ^4 

and fg are given. 
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probability density function f2(R* B;t) 
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= a Jb £ d 

P = (a - 4Es)r E0 

O b 

Q = (s^ + jEb)T E0 




Figure 4 b; The annular sector D^(R, B:ir,T;t) 
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Figure 4 c: The region 
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for 3 discrete values of 9 , namely. 



H 
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Figure 4 e: The case where is a rectangle 

0 elsewhere 





of which is "swept out". 
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CHAPTER 5 



GENERAL I^EMARKS ON TORPEDO SALVOS 
5. 1 Torpedo Lethality 

The number of T/0'§ required, on the average, to destroy a 
target, depends on two primary factors: 

(1) The accuracy of the T /O, i. e. the probability the T/O 
will hit the target, and 

(2) The destructiveness of the T/O, i. e. the probability 
that the target is destroyed if the T/O hits it. 

Here (1) is reduced to the problem of finding the I^'s , 
thence n , 

For (2) in many cases a single hit i^ not enough to destroy a 
target. Even for merchant ships a singly T/O hit is not usually 
enough to sink it, and heavy combatant s^ips are designed to 
withstand many T/O hits. To t^^at such cases exactly, one should 
determine the probabilities, say, d^ , d^ • d^ , etc. of sinking 
the ship if the ship is hit 1, 2, 3, etc. . times. Then if, for a 
given method of aiming T/O's, the prob^jkbilities of 1, 2, 3, etc. 
hits are P^ • P 2 $ » etc, the probability of sinking the ship 

would be djPj + ^ 2,^2 *^3^3 where the d*s would repre- 

sent damage coefficients. 

For this paper it is simply assumed that there is a number m 

for each type ship such that, on the average, m hits of that type will 

n 

srihk the: with probability = P(iqft,n,d)= ^ P. 
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with d s ij L sin A the spacing along the normal to T/O track* 
for the a. p. s. described in Chapter 2. 

5. 2 Time Intervals in Salvo Fire 

' f »■ ■ ■■ I * ■■ 

Recall that, in firing tiorpedoes from sL submarine, the firing 
is not a true spread in the senes that adl t)o's are launched at 
once to provide a eimultaneoue menace to t^^ enemy. T/O's are ' 
fired singly and consecutively with a prescribed time delsly to 
avoid countermining. 

Let A*y = time delay (an interval of iteveral secs. ). Consi- 
der a spread of nxm T/O's, let ms 1 . Let y' be initial 
make ready time after last inputs to TDC at time t . The mean 
T/O would go at time after time t , on an 

expected run time mean ; 

then mean =Tj +^2 *0 * ^ ^ 

The mean track To , commenced at time f + 7j be 

considered as the datum line about which the spread, or aim 
point sequence (a. p. s . )' is placed. 
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CHAPTER 6 



THE FIRING PHASE DECISION 
6, I The Cardinal Rule 

The number of aim points n in an a. p. s. may be regarded 
as depending upon the function f(t^;£) introduced in Section 4. 2, 
Hence n is a function of 'Y'. 

Therefore ( as a ojjjfdinal rule to minimize the spreading and 
thinning out of probability density, avoid protracted firings: 

FIRE AS SOON AS POSSIBLE AFTER DECISION TO 
COMMENCE FIRING PHASE. 

6. 2 The General Procedure 

Most of the notation used here has already been discussed, 
otherwise it is defined in the glossary. Four sets of procedures 
are outlined, following the dictates of the four objectives pro- 
posed in Chapter 1. In each case n ^'maximum i for which 

P. ^ p^ and d s ^ L sin A along the normal to T/O tracks 

n 

in the aim point sequence. Then P = ^ P. 

isl 

I. Fire Salvos to maximize probability of at least one hit 
(i. e. case Bl). Do not fire T/0''s at aim points for which asso- 
ciated probability of hit is less than p^ . As nearly as possible, 
all T/O's are to be expended in not more than T days; obtain 
extra hits as indicated by frequency of attack. 

(a) If 2 < n fire one T/O at each of first aim points 
(a. p. ) (cf. Section 2. 3) 
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(b) If n lookatC.' then if: 

3 

(i) 2 T/O's at each of first n a.p. and 

T/O at each of first - ^n a.p. 

(ii) n < ^ I < ^ , fire u T/O's at each of first n a, p. 

(iii) S' < n < ^ fire one T/O at each of first n a.p. 

*1 

(c) When ^ « 2]^ request permission to reduce value of p 

(d) Return to base when all T/O's are expended, but do not 
remain on station more than T days. 

II. Fire T/O's to maximize probability of m hits on targets 
requiring m hits to be sunk with probability Po (i. e. , case Al)'. 

Do not fire at aim points for which associated probability of 
hit is less than p 

As nearly as possible, all T/O's are to be expended in not 
more than T days; obtain extra hits as indicated by frequency 
of attack. 

(a) If mn sTl , fire m T/O's at each of first \) a. p. 
(Optional to fire 2- ml) T/O's at ( 1 )+ l)th a.p.). 

(b) If mn < ^ , determine '. then if: 

(i) 2^ < fire yiA. T/O's at each of first n a.p. 

(ii) mn< ^ y*-' T/O's at each of first 

n a. p. 

(iii) ^^l < mn < fire m T/O's at each of 

first n a. p. 

(c) When YJ'- ® request permission to increase 

J 

value of Po and/or to reduce value of p . 

■^o 
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(d) Return to base when all T/O's are expended, but do not 

remain on station morls thnn T. days. 

”/ 

r 

111 . Remain on station for approximately days, conducting 
as many attacks as possible. Ffre salvos to maximige probability 
of at least one hit (i. e, , case B2) 

Do not fire T/O's at aim points for which associated probabi* 
lity is less than p^ . As nearly as possible, expend all T/O's by 
end of patrol; obtain extra hits as indicated by frequency of attack. 

(a) If S fire one T/O at each of first a. p. 

(b) If n < fire T/O's at each of first n a. p. 

and one T/O at each of first ^ - / 4 .n a. p. . 

1 1 

(^) When * 2 request permission to reduce value 

of Pg- 

(d) Return to base when all T/O's are expended, but make 
this as close as possible to T^ days. 

pf. Remain on station for approximately T^ days, conduct- 
ing as many attacks as possible . Fire salvos to maximize probabi« 
lity of m hits on targets requiring m hits to sink with probability 
P (i. e., case A2). 

Do not fire T/O's at aim points for which associated probabi- 
lity of hits is less than p^ . As nearly as possible expend all 
T/O's by the end of the patrol; obtain extra hits as indicated by 
frequency of attack. 

(a) If ^ mn fire m T /O' s at each of first '\)j a. p. 

(Optional to fire 2l “ T/O's at (*i) ^ + l)_th a.p. ). 
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(b) If mn < fireyUjT/O's at each of firat n a. p. (Optional 

to fire ~ T/0‘a at (n + 1) ^ a. p. ) 

(c) When 2 j “ ^2 j permisaion to increaae value of 

P and/or reduce value of p_ . 
o o 

(d) Return to baae when all T/Q'a are expended, but make thia 
as close as possible to days. 
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Appe ndix A 

The Normal Distribution 



A-1 The One -Dimensional case. 

The normal distribution is one of the most significant mathema- 
tical models, from both the theoretical and practical standpoints, 
upon which the techniques of applied statistics are founded. 

Using as an example, say, a set of n measurements, with 

errors x^(i =1,2, ?n) , some underlying assumptions of the 

distribution are; 

(a) In .performing the set of measurements, not all the errors 
X. are equally likely to occur. 

(b) In general, large errors are less likely to occur than 
small ones. 

(c) If there is a best value m , then positive and negative 
errors about this value are equally likely to occur. 

(d) If (c) is true, then m is the arithmetic average of the 
individual measurements. 

The density of this distribution is given by 



Thus, a normal distribution is completely determined by specify- 



The function is shown in Fig. A- 1 for several values of ( 5 ~ . 
Changing m merely shifts the curves to the right or left without 
changing their shape. 




ing its mean m , and standard deviation i, e. f(x;m,(5'). 



A-1 



x-m 



By the linear transformation* y = -- , f(x) is put in stan- 
dard form, with m = o and ^ = 1 * The density function is then 

f(x;o,l) = Exp ( - 4 ) 

/ 2ir 

It is in this form that one ustially finds the integral tabulated for 
any value of x . 

A-2 The Two-Dimensional Case. 

A simple case of two independent, normally distributed varia- 
bles will be discussed here. Given f,(x;x , (T :) and f,(y;y , (T ') 

1 o X ^ o y 

such that f(x,y) = f,(x)f,(y) . Let the point (x , y ) = a 

lb o o o 



The joint density is given by: 

«*•■>'> = ii - J y {- 2 




y-a. 
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(i. e. unit mass is spread over a symmetrical mound). Fig. A-2, 
Curves formed by a vertical slice parallel to an axis are of 
the same shape as the normal distribution curve. 

The curve formed by a horizontal slice is an ellipse 

^ \ 2 ^ \ T 

f'x-a^ \ 



f(x.y) = 




= c (a constant) 



The equation represents an ellipse with center at semi- 

axes of length (c(T , c (T ) and axis parallel to (x, y) axe's. An 

X y 

ellipse of this nature is called an iso-density contour ellipse. 

The probability that (x, y) lies inside the ellipse equal-s P when 
2 2 2 

c = /\^ p , f = 2 (where c is found using the x, y values and 
P is the probability in per cent for the corresponding fractile of 
the Chi-squared distribution with 2 degrees of freedom). 
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To define a bound for the area D covered by f(x,y), any 
contour ellipse ^ maybe selected, e.g. noting that: 









90 



then c = 4. 61 



= 4.61 (f=2) 

c = 2.15 



let ^ = 90% contour ellipse with semi-axes (2. 15(T ,2. 15{T ) 

X y 

i. e. for ^ = 90% the ellipse is 4. 3(J^ by centered at a^. 




Figure A-1: Normal Frequency Curves. 

m=0; (T = 0,4,1. 0,2. 5 
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Appendix B 

Proposal for Preparation of Histograms Used in 
Salvo Firing 

A proposal for using an electronic digital computer for deter- 
mining D and F(D) suitable for use as reference diagrams. 

First, prepare' statistically, the histograms representing 
R, f(R); B, f(B); 6, f(6) where 6 = 6(C,B); and S, f(S). 

Thus for R, we might have the histogram shown in the following 
figure. I 




^1 ^3 ^^4 ^5 ^6 ^ 



F or each component the machine is instructed to pick a random 
number (by a s|lb-routine method) say 0 ^ r ^ 1 and this num- 
ber instructs the machine from which step of the cumulative graph 
to pick the value of R. e.g. , 0 r A^^ picks R^, Aj^ < r ^ 
Aj^ + A^ picks R 2 » etc. In this manner, a value of each compon- 
ent is picked and the machine computes the simple trigonometry 
involved to get the X, Y coordinates, say from the datum point 
for a given 
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These are compared to the grid values, sorted, and added to the 
proper cell as a simple count of. the occurrence 




Suggest 
100 yard 
boxes 



D and £(D) for various set of components and T^'s could 
be make up to guide the approach officer in the selection of the 
p.'s and n . 

Not only would one like histograms based on various mean 
values to use in the computation, but also on the method of 
measurement, such as range histograms at various ranges 
taken by: radar, optical, sonar, etc: Bearing: optical, hydro- 
phone, sonar; Speed; screw count; tracking, etc. 

If each histogram was made up in standard steps the compu- 
ting routine could be used' over and over again to construct D's 
and f(D) for various combinations by just changing the values 
in the cumulative histogram cells. 

Note that the computer would produce discrete mounds that 
would be as good as the histogram data from which they were 
prepared. 



B-2 



No mention of what the component error data actually looks 
like is made here, because of the non-classified nature of the 
paper. 
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